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Summary g 

This paper presents the results of an analy- 
tical study aimed at predicting the aeromechani- 
cal stability of a helicopter in ground reso- g 

nance, with the inclusion of aerodynamic forces. R 

The theoretical results are found to be in good 
agreement with the experimental results, avail- @ 

able in the literature, indicating that the u 

coupled rotor/fuselage system can be represented 

by a resonably simple mathematical model. g g 
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function 
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= profile drag coefficient of 
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= stiffness of the root springs C »C 

of the blade in flap and lag lK n ' 

respectively 
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= mass/unit length of the blade 
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= rotor radius c 


s - complex eigenvalue 
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*= blade precone, in the equations 
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= nondimensionalized quantity, 

with respect to R when involving 
length, and with respect to 
when involving frequency 


I, Introduction 

The aeromechanical instability of a helicop- 
ter, on the ground and in flight, is caused by 
coupling between the rotor and the body degrees 
of freedom. This instability is commonly denoted 
air resonance when the helicopter is in flight 
and ground resonance when the helicopter is 
on the ground. The physical phenomenon involved 
during this instability is quite complex, the 
rotor lead-lag regressing mode usually couples 
with the body pitch or roll to cause the insta- 
bility. The nature of the coupling which is both 
aerodynamic and inertial is introduced in the 
rotor due to body or support motion. Development 
of a mathematically consistent model capable of 
representing the coupled rotor/fuselage dynamic 
system is of fundamental importance for the study 
of these type of problems. The mathematical model 
should be consistent because the geometrically 
nonlinear terms associated with moderate blade 
deflections are known to have a significant role 
in rotary wing aeroelasticity^ . Thus various 
terms having the same order of magnitude must be 
retained throughout the derivation of the equa- 
tions of motion. A consistent mathematical model 
has been developed^ »3 # by the authors, to study 
the aeroelastic, structural dynamic and aero- 
mechanical effects in multi-rotor systems. 

Bousman^ has obtained excellent experimental 
data for aeromechanical stability of a hingeless 
rotor on a special gimbaled support simulating 
body pitch and roll degrees of freedom. The 
availability of this high quality experimental 
data provides an opportunity for comparing the 
results obtained from the analytical model with 
this experimental data. Bousman attributed some 
of the discrepancies found between the theoreti- 
cal results presented in his paper and experi- 
mental results to dynamic inflow. This conclu- 
sion was also examined by Johnson^ , in a recent 
study, where unsteady aerodynamic effects on the 
rotor was represented by a perturbation inflow 
model^. Johnson showed that theoretical results 
based on his model^, with dynamic inflow pro- 
vided results which showed better agreement with 
the experimental results than the results b^sed 
on a quasi-steady aerodynamic model without 
dynamic inflow. He concluded from his study 
that unsteady aerodynamic effects are repre- 
sented quite well by a dynamic inflow model. 

Using the mathematical model developed by 
the authors2»3 j it is shown that the theoretical 
results, based on the quasi- steady aerodynamic 
model, are for most cases in better agreement 
with the experimental results than the agreement 
noted by Bousman^ . The agreement with the 
experimental data is also comparable to that 
obtained by Johnson^ , except that the quasi 


steady model is incapable of predicting the 
"dynamic inflow mode" found by Johnson^, which is 
a result of the augmented state due to inflow 
dynamics . 

The good agreement between the analytical and 
experimental results indicates that the relatively 
simple analytical model is accurate for this case. 
Furthermore it also implies that only part of the 
discrepancy between theory and experiment , found 
by Bousman, may be attributed to dynamic inflow. 

II. A Brief Summary of the Experiment 

A clear description of the experimental set 
up, used for simulating the fundamental aspects of 
the aeromechanical stability of a hingeless rotor 
helicopter , is presented in Ref . 4 . The rotor 
consisted of three blades and five different con- 
figurations were tested. The different configura- 
tions represent different blade parameters char- 
acterized by the nonrotating natural frequencies 
of the blade in flap and lag, pitch-lag coupling 
and flap-lag coupling. The rotor was designed 
such that most of the blade flexibility is con- 
centrated at the root by building in root flexures. 
The rotor assembly was supported on a gimbal which 
had pitch and roll degrees of freedom. In this 
paper the analytical results obtained were com- 
pared with the experimental results, presented by 
Bousman, for rotor configurations 1 and 4, where 
the designation of these configurations is con- 
sistent with those in Bousman' s paper 4 . 

A brief description of these configurations 
is presented for the sake of completeness. Con- 
figuration 1 had different stiffnesses in flap 
and lag respectively, the corresponding nonrota- 
ting flap frequency was 3.13 Hz and that for 
lead-lag was 6.70 Hz. Configuration 4 was a 
matched stiffness case where the nonrotating flap 
frequency was 6.63 Hz and that for lead-lag was 
6.73 Hz. The airfoil cross-section of the blade 
was cambered and had a zero lift angle of attack 
equal to -1.3 degrees. Thus a substantial part of 
the experimental data was obtained for zero pitch 
setting, however, due to the presence of camber 
the rotor produces a small amount of thrust at 
this pitch setting. The rotor blades were rigid 
outboard of the flap and lag flexures which were 
located at a radial station 0.105R. There was no 
flap-pitch or pitch-lag couplings for these two 
configurations (configurations 1 and 4). Further- 
more, the blade was very stiff in torsion. In 
the case of the experiments conducted for pitch 
angles other than zero, the experimental set up 
was so designed as to introduce the changes in 
pitch angle outboard of the flexures and hence 
there was no flap-lag structural coupling for 
these cases. The structural damping in body roll 
was very small in comparison with that for body 
pitch. The body pitch and roll frequencies were 
controlled by cantilever springs on which the 
gimbal was mounted. It is stated in Ref. 4 that 
the body pitch spring was selected to provide a 
dimensionless body pitch frequency of about 0.12 
at the nominal rotor speed of 720 R.P.M. and the 
roll spring was selected to give a dimensionless 
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roll frequency of about 0.28. (The frequencies 
are nondimensionalized by dividing by rotor 
speed.) Based on these values, the dimensional 
frequencies in pitch and roll are 1.44 Hz and 
3.36 Hz respectively. It was also mentioned in 
Ref. 4 that the dimensional values of the body 
pitch and roll frequencies are about 2 Hz and 
4 Hz respectively. So the difference, noted 
between the two sets of body frequencies, raises 
a question as to what are the exact values for 
the pitch and roll frequencies. However the 
experimental results presented in Ref. 4 showed 
that over a wide range of ft (200 - 1000 R.P.M.) 
the pitch and roll frequencies are very close 
to 2 Hz and 4 Hz respectively. Hence, for the 
present analysis, the pitch and roll frequencies 
are selected to be 2 Hz and 4 Hz. The reason 
for choosing 2 Hz and 4 Hz for body frequencies 
was that at approximately 750 R.P.M. , the lead- 
lag regressing mode of the rotor was close to 
the body roll frequency causing an aeromechanical 
instability. The data used in our calculations, 
is presented in the Appendix B. 


The inflow ratio a, used in the calculation 
of the aerodynamic loads was evaluated from^ 


A = 


Oa 

16 


24 6 


I + 


eff 


aa 


where a is the solidity ratio 


( 1 ) 


a is the lift curve slope 

and 9 e ££ is the effective angle of attack of the 
blade . 

As indicated in Ref. 4, a cambered airfoil was 
used in the model rotor tested, thus 


eff G c 0 ZL 


( 2 ) 


where 9 c is the collective pitch setting of the 
blade 


is the zero lift angle of attack. 


III. Description of the Analytical Model 

The analytical model used to study this 
aeromechanical stability problem is based on the 
equations developed for a multi-rotor system 
presented in Ref. 2 and 3. Those equations 
represent the dynamics of the coupled rotor/ 
vehicle system consisting of two rotors inter- 
connected by a flexible structure. The various 
degrees of freedom considered, in deriving the 
equations, are flap, lag, torsion for each blade, 
rigid body translation and rotation of the com- 
plete vehicle and also the degrees of freedom 
representing the normal modes of vibration of 
the supporting structure. From this multi-rotor 
analytical model, only those degrees of freedom 
and the corresponding equations of motion that 
are relevant for the present study have been 
retained. The most important assumptions upon 
which the formulation is based on are: (1) the 

rotor consists of 3 or more number of blades, 

(2) the rotor is lightly loaded, (3) the rotor 
is in uniform inflow, and (4) the rotor blade 
is modelled as a rigid blade with orthogonal 
springs located at the root of the blade (Fig. 

1), where and represent the stiffness of 
the blade in flap and lag motions. 

The aerodynamic model is based on 
Greenberg's® derivation of unsteady aerodynamic 
loads on an oscillating airfoil in a pulsating 
flow. This theory is basically a modified form 
of Theodorsen's unsteady aerodynamic theory. 

By assuming the Theodorsen's lift deficiency 
function C(k) = 1 and neglecting the torsional 
motion of the blade, the aerodynamic model 
becomes a simple quasi-steady model with apparent 
mass terms. In the present calculations, only 
this quasi-steady aerodynamic model with apparent 
mass terms is used. It was found from our cal- 
culations that neglecting the apparent mass 
terms from the aerodynamic model affects the 
results only by 2 - 4%. 


The zero lift angle of attack, for the airfoil 
employed ^ (NACA 23012), was 0 £l = -1.5 degrees. 

As mentioned earlier, the equations of motion 
are nonlinear, because geometrical nonlinearities 
due to moderate deflection of the blade are 
included. Retention of the nonlinear terms is 
based upon an ordering scheme** 2 . The blade 
degrees of freedom, representing blade slopes are 
assigned an order of magnitude represented by a 
symbolic quantity e, and are denoted to be of 
order 0(e), where 0.1 < e < 0.15. The fuselage 
degrees of freedom are assumed to be of a slightly 
smaller magnitude 0(e 3 / 2 ). As indicated in Ref. 
1, this assumption is quite important for obtain- 
ing equations which are manageable from an alge- 
braic point of view. The ordering scheme consists 
of neglecting terms of order 0(e 2 ) when compared 
to order one, thus 1 + 0(e 2 ) - 1. 

The degrees of freedom considered in this 
aeromechanical stability analysis are: the fun- 

damental flap and lag modes for each blade and the 
pitch and roll degrees of freedom of the body. 

In this class of problems, it has been established 
that the collective flap and lag modes do not 
couple with the body motion and thus, these modes 
are not considered. Therefore, the total number 
of degrees of freedom governing the aeromechanical 
problem are six. These consist of: cyclic flap 

(Blc»Bi s ), cyclic lead-lag (Ci c »Gl s )» body pitch 
(0) anc body roll (<f>) . 

IV. Method of Solution and Discussion 
of Results 

The method of solution for coupled rotor/ 
fuselage problem follows essentially the procedure 
explained in Ref. 1 and 3. A brief outline of 
the procedure is given in the following few 
paragraphs . 

The equations of motion, for coupled rotor/ 
fuselage problem, are usually nonlinear 
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differential equations with periodic coeffi- 
cients. These differential equations can be 
either ordinary or partial depending on the type 
of model used for the representation of the 
rotor blade. If the blade is modelled as a 
rigid blade with root springs, the resulting 
equations will be nonlinear ordinary differential 
equations. On the other hand, if the blade is 
modelled as a flexible beam, the final equations 
will be nonlinear partial differential equations. 
In this case, the partial differential equations 
are first transformed into ordinary differential 
equations using Galerkin’s method. Thereafter the 
method of solution is the same regardless which 
of these two blade models is used. In the 
present case, because the blade is modelled as 
rigid blade with root springs (Fig. 1), the 
equations of motion are nonlinear ordinary 
differential equations with periodic coeffi- 
cents. The steps involved, in solving these 
equations to obtain the stability information, 
are as follows. 

1. Evaluation of the equilibrium position for 
the blade. 

2. Linearization of the nonlinear ordinary 
differential equations about the equilibrium 
position. (Linearized equations will have 
periodic coefficients.) 

3. Transformation of the linearized equations 
with periodic coefficients to linearized 
equations with constant coefficients, by 
applying multiblade coordinate 
transformation. 

4. Evaluation of the eigenvalues of the 
linearized system with constant coefficients 
to obtain the information about the 
stability. 

For the case of hover, the equations which 
represent the static equilibrium of the blade 
are obtained by imposing the requirement that 
all time derivatives of the blade degrees of 
freedom and the fuselage perturbations vanish in 
the equations. The resulting equations are non- 
linear algebraic equations and they are identical 
for all the blades in the rotor indicating that 
the static equilibrium is same for all blades. 
This static equilibrium position is obtained by 
solving the nonlinear algebraic equations using 
a number ical method, namely the Newton-Raphson 
technique. Then the blade degrees of freedom 
are expressed as time varying perturbations 
about the static equilibrium position, 3 q and 4 q 
for flap and lag respectively. 

3 k OIO = 3 0 + A3 k 0!0 
- C 0 +Ac k (M 

Substituting these into the nonlinear 
ordinary differential equations of motion and 
neglecting terms which contain the product of 
the perturbation terms, yields the linearized 


equations of motion. The linearized equations for 
the k-th blade will have periodic coefficients, 
since the k-th blade equations are written in the 
blade fixed rotating coordinate system. Trans- 
formation of the perturbations equations to a non- 
rotating system will result in equations with con- 
stant coefficients. This transformation is per- 
formed using the multiblade coordinate transforma- 
tion^. During this transformation, the individual 
blade degrees of freedom will transform to a new se 
of rotor degrees of freedom. In the past, these 
rotor degrees of freedom have been referred to as 
multiblade coordinates or Coleman coordinates or 
Fourier coordinates or rotor-plane coordinates. 
These coordinates are basically representative of 
the behavior of the rotor as a whole when viewed 
from a nonrotating frame. For the sake of com- 
pleteness the equations of blade equilibrium, the 
linearized perturbational blade equations (in the 
multiblade or rotor plane coordinate system) and 
the perturbational equations for the pitch and roll 
degrees of freedom are presented in Appendix A. 

Stability of the linearized system is 
determined by performing an eigen-analysis on the 
linearized constant coefficient perturbation 
equations. The eigen-values appear as complex 
pairs s = CT + iw. The complex part of the eigen 
value (w) refers the modal frequency and the real 
part (a) refers the modal damping. The mode is 
stable if C is negative and it is unstable if O 
is positive. 

For the present problem, there are six pairs 
of complex eigen-values each one representing one 
of the six degrees of freedom, namely, 3l c > 3ls> 
Clc> ^ls> 0 and The modes corresponding to 

the rotor degrees of freedom (3ic» 3ls> Clc* £ls) 
are referred to either progressing mode or 
regressing mode. The designation of progressing 
or regressing to a particular mode is based on 
the numerical value of the rotating natural fre- 
quency of the rotor. Suppose the rotating 
natural frequency, say in lead-lag, is f/rev. 

Then the two frequencies corresponding to the 
cyclic lag modes (^ic» Cls) w iH be (f+l)/ re v 
and (f-1) / rev , where f+1 is the high frequency 
lag mode and f-1 is the low frequency lag mode. 

If f is greater than 1/rev, the high frequency 
lag mode (f+1) is a progressing mode and the low 
frequency lag mode (f-1) is a regressing mode. 

On the other hand, if f is less than 1/rev, the 
high frequency lag mode is a progressing mode 
and the low frequency lag mode is also a progres- 
sing mode. These designations are also appli- 
cable for the flap modes of the rotor. A clear 
description of these is given in Ref. 6. For a 
stiff-in-plane rotor, the rotating natural 
frequency in lead-lag greater than l/rev. Hence 
the high frequency lead-lag mode is a progressing 
mode and the low frequency lead-lag mode is a 
regressing mode. For a soft inplane rotor since 
the rotating natural frequency is less than 
1/rev, both high frequency and low frequency lag 
modes are progressing modes. 
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In the present study, aimed at the aero- 
mechanical stability of a model helicopter, the 
behavior of the model is studied at various ft's 
of the rotor. Thus due to the variation in ft, 
a stiff inplane rotor at low ft's will become a 
soft inplane rotor at high ft's. In the experi- 
ment performed by Bousman^, the flexibility of 
the blade in lead-lag is such that the rotor 
becomes a soft inplane rotor beyond ft = 445 
R.P.M. Hence, for ft < 445 R.P.M., the lead-lag 
modes will have one progressing mode and one 
regressing mode and for ft > 445 R.P.M. , both the 
lag modes will be progressing modes. In Refs. 

4 and 5, even for ft > 445 R.P.M., the low- 
frequency lag mode is referred as regressing 
mode instead of progressing mode. The reason 
could be to avoid any confusion while referring 
to various modes. So, for the sake of con- 
sistency, during the discussion of our results, 
the low frequency lag mode is always referred 
as lag regressing mode. 

The results for configuration 1 are pre- 
sented in Figs. 2-7, while the results for 
configuration 4 are presented in Figs. 8-12. 

The variation of the various modal frequencies 
with ft are presented in Fig. 2, together with 
the experimental data obtained in Ref. 4. The 
progressing flap and progressing lead-lag 
frequencies increase very rapidly with ft. The 
lead-lag regressing mode frequency evaluated 
from the analytical model is in excellent 
agreement with the experimental results. The 
body pitch and roll frequencies have slightly 
higher values than the experimental results. 

The damping in pitch as a function ft is shown 
in Fig. 3. The analytical results are in 
relatively good agreement with the experimental 
data. The variation of the damping in roll as 
a function of ft is shown in Fig. 4. It is 
evident that for this case the analytical 
results yield values which are somewhat higher 
than the experimental data. The differences 
observed between our analytical results and the 
experimental points, for the frequency and damp- 
ing in body modes, could be explained as fol- 
lows. In our calculations, the numerical values 
used for the stiffness and structural damping 
in body pitch and roll modes are evaluated 
based on pitch frequency equal to 2 Hz and roll 
frequency equal to 4 Hz. As pointed out in 
Sec. II of this paper, there is some doubt 
about the correctness of the body frequencies 
(2 and 4 Hz) because in Ref. 4, there are two 
different sets of frequencies for pitch and 
roll, namely 1.44 and 3.36 Hz, and 2 and 4 Hz 
respectively. 

Figure 5 represents the variation of damp- 
ing in lead-lag regressing mode with ft. As 
indicated before, Johnson's results^ show that 
the theory with inflow dynamics shows better 
agreement with experimental data than the 
theory with quasi-steady aerodynamics. However, 
even with quasi-steady aerodynamics, the results 
of the present analysis show slightly better 
agreement than the results obtained in Ref. 5 


with inflow dynamics. It is also important to 
note that in the region, beyond 800 R.P.M., our 
results are in excellent agreement with the 
experimental data, while the theory with inflow 
dynamics predicts higher values. 

Results from the calculations performed 
indicated that the progressing and regressing 
flap modes are always stable and the damping in 
these modes increases monotonically with ft for 
configuration 1 as well as for configuration 4. 
Since these modes are always stable, the results 
are not presented in this paper. 

Changes in the damping of the lead-lag 
regressing mode as a function of the collective 
pitch setting of the blade are presented in Fig. 

6. Since Johnson^ has not presented a corres- 
ponding set of results, it was not possible to 
compare these results with an analysis based on 
the dynamic inflow model. At ft = 650 R.P.M., the 
results shown in Fig. 6a indicate that the 
theoretical analysis used by Bousman^ predicts a 
much lower value for the damping than the experi- 
mental results. The present analysis shows con- 
siderably better agreement. It should be noted 
however that for larger values of pitch setting 
the difference between the predicted results and 
the experimental results increases. This dif- 
ference could be attributed to the simple quasi- 
steady aerodynamic model used in our analysis. 

This difference however is much smaller than the 
one exhibited by Bousman’s results. Even more 
interesting are the results presented in Fig. 6b, 
corresponding to ft=900 R.P.M. For this case exper- 
imental results indicate a lead-lag regressing 
mode which is always stable, but the theoretical 
results shown by Bousman^ imply an instability 
which becomes stronger beyond a collective pitch 
setting of 2 degrees. As evident from Fig. 6b, 
the results of our analysis predict the correct 
trend and the predicted damping levels are much 
closer to the experimental results. It should be 
noted again that the agreement between the pre- 
dicted and experimental results diminishes with 
increasing collective pitch setting. An item to 
be noted in these figures (6a, 6b) is that the 
curve representing our analytical results starts 
from an angle 0 C = -1.5 degrees. Although Fig. 6 
contains an experimental data point corresponding 
to 0 C * -3 degrees, we have not computed the 
results for this pitch setting because for 0 C = 

-3 degrees, the relation between inflow ratio 
and the collective pitch of the blade (Eq. 1) 
becomes indeterminate. 

The variations in pitch damping as a func- 
tion of collective pitch setting are shown in 
Fig. 7a, and similar variations for roll damping 
are shown in Fig. 7b. As evident from 7b, the 
damping in roll is predicted quite well. However 
the damping in pitch predicted by the present 
analysis is much lower than the experimental 
results. One can only speculate on the possible 
cause for this discrepancy. One possible reason 
could be the slight nonlinearity present in the 
structural damping in pitch mentioned in Ref. 4. 
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At ft = 650 R.P.M., the lead-lag regressing mode 
frequency is close to the body pitch frequency 
(Fig. 2) and therefore the amplitudes in pitch 
could be higher. Thus nonlinearity in structural 
damping in pitch could manifest itself by 
increasing the total damping in pitch. 

The results for configuration 4 are pre- 
sented next. The variation of modal frequencies 
with ft are shown in Fig. 8. The lead-lag 
regressing mode frequency is in excellent 
agreement with our analytical predictions. The 
pitch and roll frequencies are predicted well. 
Bousman' s^ experiments showed the presence of 
a frequency of about 0.8 Hz beyond ft = 350 
R.P.M., whereas the present analysis has not 
predicted any frequency close to this value. 

Note that the regressing flap mode frequency 
is close to the pitch mode over a wide range of 
ft(400 < ft < 1000 R.P.M.). Thus it is possible 
that the pitch mode can be excited by the 
proximity of the regressing flap mode. The expla- 
nation for the presence of the 0.8Hz frequency, 
measured in the test, posed a problem since the 
theoretical results presented by Bousman^ as well 
as those obtained by Johnson^, with the quasi- 
steady aerodynamics, were incapable of predicting 
a 0 . 8Hz frequency. It is quite relevant to 
quote Bousman on this matter, Ref. 4, p. 53. 

In Bousman ’s words, ‘'However in the experimental 
case, measurements in the pitch coordinate show 
two modes of comparable damping at rotor speeds 
beyond 350 R.P.M., one mode at about 0.8 Hz 
and the other at 2 . 0 Hz". Bousman refers one 
as pitch mode (0.8 Hz) and the other as flap 
regressing mode (2.0 Hz). However, in identi- 
fying these modes Bousman states, "To call one 
mode the body mode, and the other flap regres- 
sing mode is somewhat arbitrary; the rationale 
used here is that as the blade pitch angle 
increases only one of these modes remains, and 
it is assumed to be the body pitch mode". But 
Johnson^, using the inflow dynamics model , was 
able to predict theoretically a frequency close 
to 0.8 Hz and he called it as the inflow mode 
and he identified the other frequency (2.0 Hz) 
as the pitch mode. Quoting Johnson, Ref. 5, 
p. 672, "That it is measurable (i.e., 0.8 Hz 
inflow mode) is surprising, since in fact the 
inflow variables A x and Ay do not correspond to 
real physical states of the system". He pro- 
ceeds to interpret this behavior as "the 
unsteady aerodynamics introduces behavior of 
the system, as observed in either time or 
frequency domain, that can be approximated by 
an additional oscillatory mode with low or 
moderate damping. Approximating the behavior 
by an additional mode implies then the exist- 
ence of additional states or degrees of free- 
dom of the system". Johnson also states that 
this behavior is observed only for matched 
stiffness case because "the flap regressing 
mode will be more coupled with the body motion". 
But examination of Fig. 6 in Ref. 5 (the results 
based on the theory with inflow dynamics) 
reveals that the flap regressing mode fre- 
quency is not near the body pitch frequency, 


so it is questionable whether coupling could occur 
between these two modes. In our analysis, however 
the results show that the flap regressing mode is 
close to the body pitch mode, as indicated in 
Fig. 8. Thus it appears that the interpretation 
offered by Johnson for the presence of the 0.8 Hz 
frequency mode and its designation as the inflow 
mode frequency is possible, albeit speculative. 

The variation of lead-lag regressing mode 
damping with ft is presented in Fig. 9. Again, the 
present analytical results are in closer agreement 
with the experimental results than those predicted 
by the theory with inflow dynamics. Figure 10 
and 11 show the variation of damping in roll and 
pitch modes with ft. The pitch damping is pre- 
dicted well. The roll damping is overestimated. 

The variation in damping levels of the lead- 
lag regressing mode with collective pitch angle, 
of the blade are shown in Fig. 12, for two 
different values of angular speed. It is evident 
from Fig. 12b that for the case of ft=1000 R.P.M., 
the theory used by Bousman predicts an unstable 
region beyond 0 C = 3 degrees, however the experi- 
ment indicates a stable configuration. The 
results of the present analysis are in good agree- 
ment with the experimental results. The agreement 
noted in Figs. 6 and 12, between the analytical 
results of our study and the experimental data, 
for nonzero values of collective pitch, seems to 
indicate that the discrepancy between theory and 
experiment for these cases, evident in Ref. 4, 
could be associated with the details of the math- 
ematical model and is not related to unsteady 
aerodynamic effects such as dynamic inflow. 

V. Concluding Remarks 

In this paper, the results of a theoretical 
analysis, of the aeromechanical stability of a 
hingeless rotor helicopter, are compared with the 
experimental results. Using a quasi-steady 
aerodynamic model, it was found that the results 
of the present analysis compare quite well with 
the experimental results. It is interesting to 
note that this correlation with experimental data 
appears to hold in both the region of zero collec- 
tive pitch angles considered by Johnson^ as well 
as in the nonzero range of collective pitch 
angles which was considered by Bousman^, but not 
by Johnson. Obviously the quasi steady aero- 
dynamic model is incapable of predicting the 
"dynamic inflow mode" which is caused by the 
augmented state of the system, when the dynamic 
inflow model is used. In an extension of this 
study which will include dynamic inflow, the 
physical meaning of the dynamic inflow mode will 
be reexamined. 

This study also indicates that the dis- 
crepancy between the predicted values of regres- 
sing mode lag damping and the experimental 
measurements, noted in Ref. 4, for configurations 
1 and 4, do not seem to be associated with 
dynamic inflow and are more likely to be related 
to the details of the mathematical model. 
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Furthermore the analytical model used in this 
study has the capability of simulating the 
experiment, with good accuracy, because it is 
based on the same blade model which was actually 
tested. 

Finally, it should be noted that the 
analytical model was based on an ordering scheme 
where blade slopes were assumed to be of order 
E and the fuselage rotations in pitch and roll 
were assumed to be of order £3/2, which leads 
to simplification in the equations of motion. 

The cases considered in the present study (both 
experimental and theoretical) were restricted 
so that only the linear first order terms in 
fuselage rotations were important. Thus other 
classes of problems, in which nonlinear terms 
in fuselage rotations are also exercised, have 
to be considered to determine the overall reli- 
ability of this particular ordering scheme. 
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Appendix A: Equations Used in this Study 

The equations of blade equilibrium, the 
linearized perturbational blade equation in 
multiblade coordinates, together with the per- 
turbational equations in the pitch and roll 
degrees of freedom are given below. 


Equilibrium Equations 
Flap : 

Q r- 2 . / -2 -2, . 2 C 

+ (w^ - Uip) sin G 


l 3 l 2 - 

+ Vfe} 


I 4 


+ - w 2 ) sin6 cosG + v B } 

U L r c C 4 p 


+ Vo {v h } 

7 3 77 2 _ tt 4 

+ V -3* “2 e} - v{ — 6 0 + 

E 3 _ 72 _ 

+ ^ (-A+2e0 o ) - ^ eA} = 0 


Lead-Lag 
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(A. 2) 


where 


-2 = K B 
F mft 2 R 3 
o K„ 


L mft 2 R 3 


p^abR 


; e = - ; l = 1 - e 


When there is no structural flap-lag coupling, the 
terms containing sin0 c and cos0 c must be deleted 
from the above equations as well as in the 
stability equations given below. 


ZL 


6. Johnson, W. , Helicopter Theory, Princeton where 0^ is the effective angle of attack 

University Press, Princeton, New Jersey, 

1980. 6 is the collective pitch setting of the 

blade 
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e is the zero lift angle of attack 
ZL 

Linearized Stability Equa_t ions 
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Appendix B: Rotor, Blade and Body Properties 

Rotor Geometry 

Number of blades 3 

Radius, cm 81.1 

Chord , cm 4.19 

Hinge offset, cm 8.51 

Blade airfoil NACA 23012 

Profile drag coefficient 0.0079 

Lock number 7.73 

Solidity ratio 0.0494 

Lift curve slope 2 tt 

Height of rotor hub above 

gimbal, cm 24.1 
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Blade Mass Properties 

Blade mass (to flap 


flexure) , gm 

209 


Blade mass centroid (Ref. 
flexure centerline), cm 

18.6 


Blade flap inertia (Ref. 
flexure centerline) , gm 

m 2 17.3 


Blade Frequency and Damping 

Conf. 1 

Conf .4 

Nonrotating flap freq. Hz 

3.13 

6.63 

Nonrotating lead-lag 
freq. Hz 

6.70 

6.73 

Damping in lead-lag (% 
critical) 

0.52% 

0.53% 

Body Mass Properties 



Rotary^inertia in pitch, 
gm m 

633 


Rotary inertia in roll, 
gm m^ 

183 


Body Frequency and Damping 



Pitch frequency, Hz 

2 


Roll frequency, Hz 

4 


Damping in roll (% 
critical) 

0.929% 


Damping in pitch (% 
critical) 

3.20% 





n, r.p.m. 

Figure 2. Modal Frequencies as a Function 
of ft, 0 c - 0 (Configuration 1). 
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Figure 3. Rorty Pitch Mode Damping as a 
Function of ft, 0 C = 0 (Configuration 1). 


Figure 1. Equivalent Spring Restrained 
Blade Model. 
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11, R.P.M. 


Figure 4. Body Roll Mode Damping as a 
Function of ft, 6 c = 0 (Configuration 1). 
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0 200 400 600 800 1000 

11, R.P.M. 


Figure 5. Regressing Lag Mode Damping 
as a Function of ft, 0_ = 0 (Configuration 1) . 
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-2 0 2 4 6 8 10 


0 C , deg 

(a) $1= 650 R.P.M. 



(b) 11 = 900 R.P.M. 


Figure 6. Lead Lag Regressing Mode Damping 
as a Function of 0 at (a) 650 R.P.M. and 
(b) 900 R.P.M. (Configuration 1). 








OUR ANALYTICAL RESULTS 

THEORY WITH INFLOW DYNAMICS (Ref. 5) 

A A THEORY WITHOUT INFLOW DYNAMICS (Ref. 5} 

0 EXPERIMENT (Ref. 4) 



Figure 10. Variation of Damping in Roll 
with ft, 0 = 0 (Configuration 4). 
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(a) U = 810 R.P.M. 



(b) U = 1000 R.P.M. 
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□ EXPERIMENT (Ref. 4) 


Figure 12. Variation of Damping in Lag 
Regressing Mode with 0 Q at (a) 810 R.P.M. 
and (b) 1000 R.P.M. (Configuration 4). 



Figure 11. Variation of Damping in Pitch 
with ft, 0^ = 0 (Configuration 4). 
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